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SECTION  1 
INTRODUCTION 

It  is  well  known  that  optical  emissions  resulting  from  deposi¬ 
tion  of  energy  in  the  upper  atmosphere  by  natural  and  man-made  events 
(e.g.,  aurora  and  high  altitude  nuclear  explosions)  can  serve  as  useful 
diagnostics  to  aid  in  understanding  the  detailed  physics  and  chemistry 
which  are  operative.  It  is  also  well  known  that  hot  plasmas  (temperatures 
of  several  eV  to  tens  of  eV  or  more)  which  are  formed  by  nuclear  explo¬ 
sions  may  cool  quite  effectively  by  emission  of  electromagnetic  radiation 
via  allowed  electronic  transitions.  (Hydrodynamic  expansion,  of  course, 
produces  a  great  deal  of  cooling  also.)  However,  once  the  temperature, 
or  more  precisely  the  electron  temperature  of  such  a  plasma,  has  dropped 
to  about  an  electron  volt  or  so,  electron  impact  excitation  of  atmospheric 
constituents  effectively  is  limited  to  low  lying  electronic  states  (and 
vibrational  states  in  molecules).  For  the  dominant  atmospheric  species, 
these  states  are  metastable,  so  "rapid"  radiation  cooling  is  shut  off.  At 
this  point,  interest  in  radiation  as  a  plasma  cooling  mechanism  tradition¬ 
ally  has  diminished  quite  rapidly,  although  spectroscopic  interest 
remained  strong. 

The  purpose  of  this  paper  is  to  illustrate  that  warm  plasmas 
(-1  eV  temperature)  at  high  altitudes  (where  interactions  with  the 
neutral  atmosphere  are  minimal)  can  rid  themselves  over  a  period  of  hours 
of  considerable  excess  energy  by  forbidden  radiation.  Such  slow  but 
steady  cooling  turns  out  to  be  significant  in  the  case  of  nuclear  explo¬ 
sion  produced  plasma.  In  addition,  the  partition  of  energy  between 
thermal  and  excited  metastable  s’  is  .jlored. 


Spectroscopic  and  other  studies  of  emissions  from  the  upper 
atmosphere  have  aided  understanding  of  a  variety  of  physical  phenomena. 
For  example,  thermospheric  temperatures  and  neutral  winds  have  been 
inferred  from  Doppler  measurements  of  01  6300  A  radiations  1_3,  precipi¬ 
tating  electron  energetics  in  the  auroral  region  have  been  deduced  from 
airglows  at  ionospheric  altitudes4,  and  F-region  plasma  depletions5  and 
propagating  waves6*7  have  been  observed  by  01  6300  A  emissions.  A 
recent  review  of  optical  processes  as  they  apply  to  the  F-region,  particu¬ 
larly  at  high  latitudes,  is  available8. 

In  the  realm  of  nuclear  weapons  effects,  optical  phenomena 
induced  by  high  altitude  nuclear  detonations  have  produced  both  dazzling 
displays  and  scientific  diagnostic  information.  In  1962,  a  1.4  meraton 
nuclear  device  (Starfish)  was  detonated  at  400  km  altitude  near  Johnston 
Island  in  the  Pacific.  As  this  was  an  announced  test  with  both  scientific 
and  military  objectives,  numerous  groups  stationed  equipment  at  a  variety 
of  locations  across  the  Pacific.  Spectroscopic  observations9*10  of 
optical  rotation  both  immediately  following  and  several  minutes  after 
detonation  confirmed  the  presence  of  allowed  and  forbidden  radiation  from 
atmospheric  species  in  the  immediate  vicinity  of  the  burst  point  and  in 
regions  of  optical  activity  ("artifical"  aurora,  etc.)  well  removed  from 
the  burst  point.  A  review  article  by  Hoerlin11  contains  color  photos  of 
Starfish  which  illustrate  some  of  the  striking  optical  features  which 
produced  the  spectroscopic  records  of  references  9  and  10. 

This  report  attempts  to  provide  reasonable  estimates  of  radia¬ 
tion  cooling  processes  which  operate  over  extended  periods  (i.e.,  hours) 
in  such  high  altitude  plasmas.  Forbidden  radiations  arising  from  radia¬ 
tive  decay  of  low-lying  metastable  states  of  0+,  0,  N+,  N,  and  02  have 
been  computed  under  the  assumption  that  the  low-lying  states  have  been 
populated  by  electron  impact.  Quenching  of  excited  states  by  electron 
collisions  has  been  included,  but  heavy  particle  quenching  has  been 
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ignored.  Time  dependent  rate  equations  for  the  ground  and  excited  states 
of  the  five  emitting  species  have  been  numerically  integrated  to  obtain 
self-consistent  solutions  free  from  assumptions  of  equilibrium  excited 
state  populations.  Section  2  of  the  report  describes  the  radiating  sys¬ 
tems  of  interest  and  presents  the  relevant  equations.  Section  3  contains 
the  results  of  detailed  calculations  which  span  radiation  periods  of 
several  hours,  and  Section  4  presents  a  discussion  of  results  and  some 
conclusions. 


SECTION  2 


RADIATION  COOLING  EQUATIONS 

Each  of  the  atmospheric  species  of  interest  here  (0+,  0,  N+,  N, 
and  02 )  possesses  a  series  of  low-lying  metastable  states  which  can  be 
populated  by  electron  impact.  These  states  are  displayed  schematically  in 
Figure  1.  In  order  to  simplify  the  calculations  a  bit,  each  set  of 
closely  spaced  (doublet  or  triplet)  states,  shown  in  the  figure,  has  been 
lumped  together  into  a  composite  stated  Thus,  for  the  calculations 
described  here,  each  of  the  species  is  modelled  as  a  system  with  two 
excited  states  above  the  ground  state. 

These  excited  states  can  be  populated  (and  de-popul ated)  by 
collisions  with  electrons,  processes  for  which  the  cross-sections  seem  to 
be  reasonably  well  known.  (Collisions  with  other  particles  also  can  alter 
excited  state  populations.  More  about  this  later.)  For  the  present  work, 
it  has  been  assumed  that  radiation  emitted  from  a  plasma  volume  element  is 
not  reabsorbed  and  that  the  total  radiation  field  is  sufficiently  weak 
that  stimulated  emission  is  insignificant.  Under  these  conditions,  the 
rate  equations  describing  the  excited  state  populations  of  a  generic  three 
level  system  are 

=  -koinoV  ko2noV  Al°ni  +  k i o n i ne+  k2o"2"e+  A20n2  (1) 


Present  use  of  the  word  "state"  is  somewhat  unprecise.  In  precise 
spectroscopic  terminology,  "levels"  are  being  combined  into  "terms." 
See  page  122  of  E.  U.  Condon  and  G.  H.  Shortley,  The  Theory  of  Atomic 
Spectra,  Cambridge  University  Press  (1964). 


-jj:1  =  k0in0ne-  k^n^-  A^nx  +  k2in2ne-  k12n1ne+  A21n2  (2) 

— 2  =  k  n  n  -  k„„n  n  -  A„  n  +  k,  n  n  -  k„,n  n  -  A„ ,  n  (3) 

dt  02  o  e  20  2  e  20  2  12  l  e  21  2  e  21  2  '  ' 

In  this  notation,  kjj  is  the  collisional  excitation  (or  deexcitation) 
rate  coefficient  from  state  i  to  state  j,  A-f j  is  the  Einstein  coeffi¬ 
cient  for  spontaneous  emission,  nj  is  the  number  density  of  particles  in 
state  i,  and  ne  is  the  electron  density.  The  ground  state,  first 
excited  state,  and  second  excited  state  are  denoted  by  subscripts  0,  1, 
and  2,  respectively.  Because  photon  absorption  and  stimulated  emission 
are  assumed  insignificant  in  the  present  work,  terms  containing  the 
Einstein  B  coefficient  are  absent  from  these  rate  equations. 

The  electron  impact  excitation  and  deexcitation  rate  coeffi¬ 
cients  used  for  0+,  0,  N+,  and  N  were  from  Ali12  and  Ali,  et  al13.  These 
coefficients  were  obtained  by  integrating  cross-section  data  with  an  elec¬ 
tron  Boltzmann  distribution.  (Other  rate  coefficients  or  cross-sections 
that  one  might  want  to  consider  are  available. )14-16  Note  that  the 
excitation  rate  coefficient,  k i j ,  is  related  by  detailed  balance  to  the 
deexcitation  rate  coefficient,  k j \ ,  by 


-4Efj/kT 


9ik1j  ■  9jkj1e 


where  the  g's  are  statistical  weights17,  a E i j  is  the  i-j  energy  level 
spacing,  k  is  thd  Boltzmann  constant,  and  T  is  the  electron  temperature. 

Excitation  rate  coefficients  for  transitions  from  the  ground 
state  (X3I~)  to  the  first  excited  (a1^  )  and  second  excited  (bl£+)  states 
of  02  were  computed  by  integrating  excitation  cross-sections  multiplied  by 


electron  velocity  over  a  Boltzmann  energy  distribution: 


k  =  -  2it 
excit  m  UkT)3/2 


f  E  o(E) 
o 


(5) 


where  m  is  the  electron  mass.  The  cross  section  for  ground  to  first 
excited  state  transitions  was  assumed  to  be 


The  cross-section  for  ground  to  second  excited  state  transitions  was 
assumed  to  be 


a 


=  3.09  x  10-17 

+  4.67xl0-^(l_!)2'1398  [l-(^)°'4]3 


(cm2) 


(7) 


In  Equations  6  and  7,  the  energy,  E,  has  units  of  electron  volts.  The 
square  bracketed  terms  are  to  be  interpreted  as 

1  -  |o  ,  E  >  E0 


These  cross  sections  are  due  to  Archer18  and  result  from  his  analytic  fit 
to  cross-section  data  of  Trajmar,.  et  al.19  Deexcitation  rate  coefficients 
for  these  transitions  were  computed  with  the  aid  of  Equation  4.  Unfortu¬ 
nately,  no  suitable  cross-sections  or  excitation  rate  coefficients  for 
first  to  second  excited  state  transitions  appear  to  be  available.  There¬ 
fore,  the  rate  coefficients  for  excitation  and  deexcitation  through  this 
transition  are  assumed  to  be  zero. 
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The  probabilities  (Einstein  A  coefficients)  used  for  radiative 
decay  of  the  metastable  states  of  0+,  0,  N+,  and  N  are  due  to  Wiese,  et 
al17.  The  reader's  attention  is  called  to  Appendix  B  of  that  publication 
wherein  rules  are  given  for  combining  transition  probabilities  for  members 
of  a  multiplet  into  a  composite  transition  probability.  Transition  proba¬ 
bilities  for  02  (0-0  band)  are  due  to  Bader,  et  al20  (a1^  +  X3I~  trans¬ 
ition),  and  from  the  life  time  data  assembled  by  McGowan,  et  al2  (b1^"  ♦ 
X3r  transition). 

Transition  probabilities,  statistical  weights,  etc.,  for  the  five 
systems  of  interest  are  summarized  in  Table  1.  Electron  impact  excitation 
and  deexcitation  rate  coefficients  are  summarized  in  the  Appendix. 

Once  the  excited  state  populations  have  been  determined,  the 
instantaneous  radiation  rates  can  be  computed  as  the  products  of  the  state 
densities  and  the  transition  probabilities.  Thus,  the  rate  at  which 
energy  is  lost  from  the  plasma  is  obtained. 

If  electron  thermal  energy  were  the  only  energy  reservoir  to  be 
tapped,  then  our  model  for  radiation  cooling  would  be  reasonably 
complete.  However,  at  the  plasma  densities  of  interest  here,  ion-electron 
collisions  are  sufficiently  frequent  to  insure  some  measure  of  thermal 
energy  transfer,  with  the  result  that  radiation  cooling  indirectly  removes 
at  least  some  of  the  ion  thermal  energy.  In  order  to  compute  this 
process,  collisional  heat  transfer  between  ions  and  electrons  was  assumed 
to  obey 

It  <iNekV  ’  Wle<kTr  KTC>  (9) 

where 


h.  =  3.6  x  10"9  (kT  )'1  “5 


cm3 /sec. 


(10) 


'•b  V 


Here  the  subscripts  refer  to  ions  and  electrons;  the  N's  are  number  densi¬ 
ties  (cm-3);  kTi  and  kTe  are  in  eV.  The  heat  transfer  coefficient, 
hie,  due  t0  Ki lb.  Stagat,  and  Stoeckly22. 

In  principle,  ion-neutral,  electron-neutral,  and  electron- 
vibrational  mode  (of  neutral  molecules)  heat  transfers  also  should  be 
computed.  Processes  involving  ion-neutral  and  electron-neutral  collisions 
can  be  evaluated  with  the  aid  of  heat  transfer  coefficients  from  reference 
22.  For  the  purpose  of  estimating  relative  rates,  a  700°K  Jacchia  atmos¬ 
phere23  has  been  used.  For  altitudes  above  200  km,  electron-neutral  heat 
transfer  is  negligible.  In  the  300  km  altitude  regime,  and  at  the  densit¬ 
ies  considered  here,  ion-neutral  transfer  roughly  is  on  a  par  with  ion- 
electron  transfer.  With  increasing  altitude,  heat  transfer  involving 
neutrals  diminishes  as  the  neutral  density  decreases  exponentially. 
Energy  transfer  between  electrons  and  the  vibrational  modes  of  nitrogen 
molecules  can  be  computed  also.  If  we  assume  the  vibrational  state 
distribution  can  be  characterized  by  a  temperature  Tv,  then  Equation  9 
(with  Ti  replaced  by  Tv)  can  be  used  provided  we  adjust  the  rate  coef¬ 
ficient.  Thus,  instead  of  hje,  we  have  used  an  electron  -  N2  vibra¬ 
tional  modes  heat  transfer  coefficient  which  is  based  on  vibrational  exci¬ 
tation  cross  sections  due  to  Schulz.24  This  coefficient  is 


hey  =  2.0  x  IQ"10  (kTe  +  0.3)  exp(-0.3/kTe) 


+  8.3  x  io”9  (kT  +  2.3)  exp(-x)  (cm3/sec) 


where 


— 2.3/kTg  ,  To  ^  Tw 

x  =  (12) 

-2.0/kTv  ,  Te  «  Tv  v  ' 

with  kTe  and  kTv  in  electron  volts.  One  finds  that  below  about  350  km 
(again  using  the  neutral  atmosphere  of  reference  13)  electron  -  N2  colli¬ 
sions  can  be  quite  significant. 
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For  the  present  work,  the  effects  of  ion-electron  heat  transfer 
collisions  as  prescribed  by  Equations  9  and  10  have  been  included.  The 
other  collisional  energy  transfer  processes  have  been  ignored,  in  order  to 
avoid  undue  complexity.  It  must  be  recognized  that  a  complete  calculation 
of  a  high  altitude  plasma  must  include  hydrodynamics,  chemistry,  etc., 
processes  which  we  are  not  prepared  to  include  in  the  present  work.  We 
leave  that  to  full  blown  nuclear  burst  simulation  which  should  include 
radiation  cooling  equations  and  all  the  relevant  energy  transfer 
processes. 

Electron  impact  collisions  are  not  the  only  means  by  which  meta¬ 
stable  state  populations  may  be  changed.  Excited  state  populations  also 
can  be  altered  (quenched)  by  collisions  with  heavy  particles,  but  results 
of  such  collisions  may  depend  strongly  on  the  initial  and  final  states  of 
the  heavy  particle.  Given  the  large  variety  of  potential  target  particle 
state  and  quenching  heavy  particle  state  combinations  which  are  possible, 
one  should  not  be  surprised  that  cross-section  or  rate  coefficient  data 
for  heavy  particle  quenching  are  rather  spotty  and  tend  to  be  fairly 
uncertain.  If  one  uses  the  limited  set  of  quenching  rate  coefficients 
suggested  by  Gerard8  or  A 1  i ,  et  al,13  and  a  nominal  neutral  atmosphere,23 
and  if  one  looks  for  the  altitude  above  which  radiative  decay  of  a  state 
exceeds  heavy  particle  quenching,  then  one  finds  that  the  present  calcula¬ 
tions  generally  apply  above  200-300  km  altitude,  the  precise  altitude 
depending  on  the  particular  excited  species  -  quenchant  combination. 

For  the  present,  heavy  particle  quenching  has  been  ignored 
also.  However,  the  lack  of  such  quenching  is  consistent  with  the  decision 
to  ignore  heat  transfer  processes  which  would  tap  the  neutral  thermal  and 
N2  vibrational  energy  pools.  Quenching  collisions  involving  neutrals 
would  tend  to  return  energy  from  the  metastable  states  to  neutral  thermal 
where,  by  heat  transfer,  it  could  find  its  way  back  to  the  electrons. 


SECTION  3 
RESULTS 

This  section  contains  results  from  a  systematic  set  of  calcula¬ 
tions  which  was  designed  to  illustrate  the  behavior  of  the  various  three 
level  metastable  systems  under  circumstances  representative  of  conditions 
in  the  ionosphere/magnetosphere  following  a  high  altitude  nuclear  burst. 
The  results  are  divided  into  two  categories;  the  first  illustrates  single 
species  systems  collisionally  excited  by  electrons  at  densities  of  104 , 
106,  and  108  cm-3.  The  second  category  represents  calculations  of  multi¬ 
species  systems  based  on  conditions  found  in  selected,  detailed  calcula¬ 
tions  of  high  altitude  nuclear  bursts. 

Computational  results  are  presented  in  four  types  of  plots.  The 
first  type  of  plot  shows  ground  and  excited  state  population  densities  as 
a  function  of  time.  In  most  cases,  the  excited  state  densities  remain 
less  than  a  few  percent  of  the  total  species  density,  so  the  ground  state 
density  curve  usually  appears  as  a  horizontal  line  across  the  upper 
portion  of  this  type  of  graph.  The  discerning  eye,  however,  will  find 
that  the  ground  state  curve  is  not  simply  a  constant. 

The  second  type  of  graph  displays  radiation  cooling  rates  as  a 
function  of  time.  The  third  type  of  graph  shows  electron  and  ion 
temperatures  as  a  function  of  time.  At  low  ion  and  electron  densities 
(104  cm'3),  one  frequently  finds  that  col  1 isional  thermal  energy  transfer 
is  slow  compared  to  coll  isional  excitation  and  radiation,  so  the  ion  and 
electron  temperatures  can  depart  significantly  from  each  other. 


The  relative  partition  between  excited  state  and  thermal  ener¬ 
gies  as  a  function  of  time  is  contained  in  the  fourth  type  of  plot.  This 
class  of  graph  displays  the  ratio  of  energy  tied  up  in  excited  states  to 
available  energy.  The  available  energy  is  defined  as  the  sum  of  excited 
state  and  electron  thermal  energies,  plus  ion  thermal  energy  provided  that 
thermal  energy  transfer  (Equation  9)  is  included  in  the  calculation. 
(Equation  9  was  included  in  all  but  two  cases  which  are  noted.) 

Because  the  objective  of  this  report  was  to  emphasize  energy 
partition  and  radiation  cooling  over  extended  periods,  most  calculations 
covered  an  interval  of  30,000  seconds  (8  1/3  hours).  In  order  to  insure 
legibility  of  plots  which  might  span  many  decades,  most  graphs  have  been 
limited  to  no  more  than  six  decades.  As  a  result,  some  curves  (e.g., 
certain  excited  state  densities)  of  particularly  small  magnitude  rapidly 
disappear  or  don't  appear  at  all.  This  plotting  effect  usually  occurs,  if 
at  all,  in  low  electron  density  calculations. 

The  calculations  presented  here  are  summarized  in  Table  2.  In 
all  cases,  initial  excited  state  densities  were  assumed  to  be  zero.  As 
Table  2  indicates,  two  categories  of  calculations  were  done.  The  first 
category  is  composed  of  a  single  species  systems  with  the  species  and 
electron  densities  set  equal.  Electron  densities  were  set  to  104,  106,  or 
108  cm*3  to  cover  the  range  of  interest  in  late-time  high  altitude  nuclear 
burst  problems;  electron  and  ion  temperatures  were  initialized  to  one 
electron  volt. 

The  second  category  of  cases  is  composed  of  multi-species  sys¬ 
tems  which  are  representative  of  a  range  of  conditions  found  in  detailed 
two-fluid  numerical  simulations  of  high  altitude  nuclear  bursts.  Species 
densities,  electron  density,  and  ion  electron  temperatures  were  as  indica¬ 
ted  in  Table  2. 
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It's  worthwhile,  perhaps,  to  note  that  excited  state  densities 
and  radiation  rates  are  not  trivially  related  to  the  electron  density. 
If,  for  example,  we  assume  the  quasi-steady  state,  then  Equations  2  and  3 
yield  the  ratio  of  first  to  second  excited  state  densities  for  our  generic 
three  level  system: 

ni  km  +  ko  n  +  Aon  +  Aoi^  +  kno  /koi  n  +  Ao  i 


n2  *02  ^|*io  +  *12J  ne  +  +  *01  ^12  ne 

Even  in  the  quasi-steady  state,  the  A's  are  not  negligible  in  general,  so 
the  population  ratio  and  radiation  rate  ratios  depend  in  a  non-trivial 
manner  on  the  electron  density.  This  is,  of  course,  already  well  known. 


The  first  series  of  calculations  presented  here  deals  with  fully 
ionized  oxygen  plasmas.  In  order  to  dispense  at  the  outset  with  questions 
about  the  significance  of  collisional  thermal  energy  transfer  between  ions 
and  electrons,  cases  without  and  with  the  transfer  process  are  presented 
in  Figures  2,  3  (without)  and  Figures  4-6  (with).  Thus,  if  one  compares 
Figures  2  and  4  (electron  density  of  104  cm-*3),  it  is  seen  that  the 
0  (ZD  )  densities  are  virtually  identical  during  the  first  10,000  seconds 
and  then  gradually  deviate  until  the  curve  in  the  absence  of  energy 
transfer  (Figure  2)  is  about  20%  below  the  curve  with  the  transfer  (Figure 
4).  A  comparison  of  radiation  rates  in  these  two  figures  shows  the  same 
trend,  e.g.,  at  30,000  seconds,  the  2D°-*-4S°  transition  is  cooling  about 
20%  slower  in  Figure  2  than  in  Figure  4.  Both  of  these  observations  are 
readily  understood  by  comparing  the  electron  temperature  curves  of  Figures 
2  and  4.  During  the  first  few  thousand  seconds  both  curves  compare  favor¬ 
ably,  and  the  ion  temperature  curve  in  Figure  4  indicates  very  little 
energy  transfer  from  ions  to  electrons.  At  30,000  seconds,  however,  the 
electron  temperatures  differ  by  about  10%  and  the  ion  temperature  (Figure 
4)  has  dropped  by  about  a  tenth  of  an  eV. 
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Figure  2.  Radiation  properties  of  0+  at  an  electron  density  of  104  err 
Ion-electron  heat  transfer  (Equation  9)  is  absent  from  this 
calculation. 


*1— 1— I— I— I— J_ t— I 

0 


1  1  i 


1 0000  '  20000 

TIME  (seconds) 


Figure  3  (Continued) 
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Finally,  the  curves  denoting  the  fraction  of  available  energy 
which  is  tied  up  in  excited  states  may  be  compared.  Both  curves  peak  at 
about  10,000  seconds,  but  the  curve  in  Figure  4  peaks  at  a  lower  value 
(even  though  the  electron  temperature  and  excited  state  densities  are 
higher  than  found  in  Figure  2)  because  the  available  energy  has  been 
doubled  by  coupling  to  the  ions.  The  most  significant  point  here  is  that 
in  each  case  the  fraction  of  available  energy  tied  up  in  excited  states  is 
a  few  percent. 

With  the  aid  of  Figures  3  and  6  one  can  make  a  parallel  set  of 
comparisons  at  the  other  electron  density  extreme  considered  here 
(108  air3).  As  one  would  expect,  the  case  with  ion-electron  energy  trans¬ 
fer  (Figure  6)  maintains  a  slightly  (-30%)  higher  electron  temperature 
than  the  case  without  (Figure  4),  so  exicted  state  densities  and  radiation 
rates  are  also  somewhat  larger.  Unlike  the  preceding  pair  of  cases,  the 
ion  and  electron  temperatures  are  now  locked  together  by  collisional  coup¬ 
ling  which  transfers  energy  to  the  electrons  faster  than  it  can  be  pumped 
into  excited  states  and  radiated  away.  In  Figure  3,  the  energy  fraction 
in  excited  states  very  rapidly  peaks  at  about  10%  while  in  Figure  6  the 
corresponding  peak  is  at  about  8%. 

These  comparisons  (and  similar  comparisons,  not  presented,  for 
the  other  emitting  species)  lead  us  to  observe  that  for  the  range  of  elec¬ 
tron  and  species  densities  considered  (where  the  electron  and  species 
densities  are  equal)  failure  to  couple  ion  and  electron  thermal  energy 
reservoirs  results  at  most  in  a  modest  alteration  of  the  excited  state 
densities,  radiation  rates,  and  energy  partition.  Now  that  the  magnitude 
of  such  coupling  has  been  established,  we  proceed  to  the  remaining  calcu¬ 
lations,  all  of  which  included  the  coupling. 

One  can  now  consider  Figures  4,  5,  and  6  as  a  set  of  calcula¬ 
tions  which  illustrate  the  electron  density  dependence  of  radiation  cool¬ 
ing  by  the  0+  system.  This  series  of  figures  shows  that  at  an  electron 


density  of  104  cm-3  it  takes  roughly  10,000  seconds  to  reach  a  quasi- 
steady  distribution  of  excited  states,  radiation  rates,  etc.  At  106  cm-3 
density,  the  time  required  to  reach  the  quasi-steady  state  is  about  200 
seconds;  at  108  cm"3  density,  the  quasi-steady  state  is  attained  "instant¬ 
ly"  as  compared  to  the  time  scale  of  the  plots.  (Note  that  the  ion-elec¬ 
tron  energy  transfer  rate  can  alter  the  scaling  of  this  time  constant 
because  the  electron  temperature  can  remain  higher  for  a  longer  time.)  As 
one  would  expect,  the  high  ion-electron  density  case  radiates  much  more 
rapidly  than  the  lower  density  cases.  In  addition  note  that  at  the  104 
cm"3  density,  the  dominant  transition  is  2D°+4S°.  At  higher  densities, 
the  2PW  and  2P°+4S°  transitions  dominate.  This  switch  results  from 
the  shift  in  relative  populations  of  the  exicted  states,  as  expressed  in 
Equation  13,  which  occurs  as  the  electron  density  is  altered.  It  is  the 
interplay  of  excitation  rates  and  transition  probabilities  which  is 
primarily  responsible  for  this  effect. 

As  noted  earlier,  the  curves  denoting  the  fraction  of  available 
energy  tied  up  in  excited  states  all  remain  in  the  range  of  a  few  percent 
or  less.  Thus,  from  the  ion  and  electron  temperature  curves  we  conclude 
that  the  time  required  to  radiate  one-half  the  energy  is  about  6,000 
seconds  at  108  cm-3,  about  30,000  seconds  at  106  cm-3,  and  (very)  roughly 
2  x  lo5  seconds  at  104  cm"3 . 

The  excitation  and  radiation  characteristics  of  nitrogen  ions  at 
the  electron  densities  of  interest  are  displayed  in  Figures  7,  8,  and  9. 
Upon  inspection  of  these  figures,  one  notes  that  the  excited  state  densi¬ 
ties  and  radiation  rates  attain  a  quasi-steady  state  virtually  instantly 
as  compared  to  the  30,000  second  time  scale  of  the  plots.  Only  the  case 
at  104  cm-3  electron  density  (Figure  7)  shows  a  perceptible  period  (a  few 
hundred  seconds)  when  departures  from  a  quasi-steady  distribution  of 
states  are  evident. 


FRACTION  OF  ENERGY  IN  EXCITED  STATES 


It's  worth  noting  the  variation  with  electron  density  of  the 
ratio  of  radiation  rates  for  the  1D+3P  and  1S-*-1D  transitions.  For  the  104 
cm-3  and  106  cm-3  electron  density  cases,  this  ratio  is  always  greater 
than  one.  At  108  cm-3  electron  density  and  while  the  electron  temperature 
is  above  -0.53  eV  (during  the  first  -1500  seconds),  the  ratio  is  less 
than  one.  This  behavior  can  be  understood  in  terms  of  Equation  13.  Based 
on  trasition  probabilities  and  photon  energies,  one  can  compute  that  the 
ratio  of  *D  to  3S  state  densities  needs  to  be  -  233  or  less  if  the 
radiation  rate  for  the  transition  is  to  dominate.  At  low  electron 
density,  electron  empact  excitation  is  insufficient  to  generate  a  ratio 
this  small.  One  observes  in  Figure  9  (early  time)  that  the  *D  to  3S  state 
density  ratio  is  about  100  (electron  density  is  108  cm-3),  but  in  Figure  8 
the  ratio  has  jumped  to  roughly  1000  (electron  density  is  106  cm-3). 

As  in  the  case  for  radiation  from  oxygen  ions,  the  nitrogen  ion 
calculations  indicate  that  energy  associated  with  excited  states  is  at 
most  a  few  percent  of  the  available  energy.  Thus,  the  temperature  curves 
allow  us  to  readily  calculate  that  a  bit  more  than  2000  seconds  is 
required  to  radiate  one-half  the  available  energy  at  108  cm-3;  about  5500 
seconds  are  required  at  106  cm-3.  One  also  notes  that  at  the  104  cm-3 
electron  density,  excitation  and  radiation  is  able  to  remove  energy  from 
the  electrons  faster  than  it  can  be  replaced  by  transfer  from  the  ions. 
As  the  higher  electron  densities,  this  is  not  the  case. 

The  radiating  characteristics  of  oxygen  atoms  over  the  electron 
density  range  of  interest  are  presented  in  Figures  10,  11,  and  12.  The 
reader  is  reminded  that  the  radiating  species  density  has  been  set  equal 
to  the  electron  density  (and  also  to  a  non-radiating  ion  density).  This 
rather  arbitrary  choice  influences  quite  strongly  the  magnitudes  of 
various  results.  The  particular  choice  of  equal  densities  was  made 
because  it  provides  a  reasonable  balance  between  the  heat  capacity  of  the 
ion-electron  gas  and  energy  stored  in  excited  states.  From  these  results 


the  reader  can  extrapolate,  in  a  limited  sense,  to  other  density  combina¬ 
tions.  If  the  oxygen  atom  density  had  been  substantially  larger  than  the 
electron  and  ion  densities,  then  most  of  the  available  energy  would  have 
gone  mainly  into  the  lowest  excited  state  at  the  same  time  the  electron 
and  ion  temperatures  dropped  to  low  values.  The  resultant  curves  would  be 
relatively  flat  and  uninteresting.  (Clearly,  the  best  approach  is  to 
compute  for  the  precise  conditions  of  interest  in  a  particular  problem.) 

The  first  point  to  note  in  this  series  of  three  figures  is  that 
the  quasi-steady  state  is  obtained,  even  at  the  lowest  electron  density, 
rather  quickly  as  compared  to  the  time  scale  of  the  plots.  However,  at 
104  cm-3  electron  density  (Figure  10),  the  radiation  cooling  and  temper¬ 
ature  curves  are  quite  flat  as  a  result  of  relatively  low  excitation 
rates.  One  sees  that  less  than  a  percent  of  the  available  energy  is  tied 
up  in  excited  states,  so  it  readily  can  be  inferred  from  the  temperature 
curves  that  only  about  6  percent  of  the  original  energy  has  been  radiated 
over  the  30,000  second  period. 


The  calculation  at  106  cm-3  electron  density  shows  more  inter¬ 
esting  temporal  behavior  and  achieves  substantial  cooling  over  30,000 
seconds.  Once  again,  the  fraction  of  available  energy  which  is  tied  up  in 
excited  states  is  quite  small  (<  2%),  so  we  readily  estimate  from  the 
temperature  curves  that  one-half  the  initial  energy  was  radiated  by  about 
11,000  seconds.  Note  that  radiation  from  the  long  lived  lD+3P  transition 
is  the  dominant  loss  channel,  a  result  which  is  not  too  surprising  in 
light  of  the  fact  that  the  XD  state  has  at  least  three  orders  of  magnitude 
greater  population  density  than  does  the  *S  state.  One  might  expect  the 
relatively  short  lived  transition  to  dominate,  but  the  ratio  of 

transition  probabilities  to  lD+3P)  is  only  about  200  and  can't 

compensate  for  the  (at  least)  three  orders  of  magnitude  difference  in 
populations  (Refer  to  Table  1). 
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Inspection  of  results  corresponding  to  an  electron  density  of 
108  cnr3  (Figure  12)  reveals  that  at  early  times  the  state  density  is 
about  two  orders  of  magnitude  smaller  than  the  state  density.  In  view 
of  the  ratio  of  transition  probabilities,  it's  not  surprising  to  find  that 
initially  the  1S-*-1D  transition  dominates.  Note  that  the  rapid  drop  in 
temperature  reduces  the  *S  state  density  and  the  1S-*-1D  radiation  rate  so 
that  after  several  hundred  seconds  the  1D+3P  transition  dominates.  This 
shift  from  predominantly  green  to  predominantly  red  radiation  is  consis¬ 
tent  with  observations10  of  the  Starfish  high  altitude  nuclear  explosion. 
(One  should  not  expect  these  simplified  calculations  to  accurately  repro¬ 
duce  the  temporal  behavior  observed  in  the  actual  nuclear  burst  event 
because  too  many  relevant  aspects  of  the  magnetohydrodynamics,  air  chemis¬ 
try,  etc.  have  been  omitted  here.)  Finally,  we  note  that  in  this  calcula¬ 
tion  one-half  the  initial  energy  was  radiated  during  the  first  2500 
seconds. 

We  now  turn  our  attention  to  radiation  from  atomic  nitrogen; 
Figures  13,  14,  and  15  contain  computational  results  at  the  selected  elec¬ 
tron  densities.  The  reader  is  reminded  that  as  in  the  atomic  oxygen 
calculations,  nitrogen  atom,  electron,  and  non-radiating  ion  densities  are 
equal.  Figure  13  presents  the  results  for  an  electron  density  of  104 
cm-3.  Note  that  the  excited  state  population  densities  never  really  reach 
a  quasi-steady  distribution,  i.e.,  the  excited  state  densities  continue  to 
increase  as  the  electron  temperature  falls  during  the  entire  30,000  second 
interval.  Also  note  that  at  zero  time  the  2P°  state  density  "instantly" 
jumps  to  a  value  near  10"1  cm-3,  and  the  radiation  rates  for  transitions 
which  originate  from  this  level  similarly  jump  to  finite  values.  A  plot 
with  a  greatly  expanded  time  scale  shows  that  the  2P°  state  density  rises 
with  an  e-folding  time  of  about  10  seconds,  so  the  rise  appears  instant¬ 
aneous  on  a  30,000  second  time  scale.  Finally,  the  fraction  of  available 
energy  tied  up  in  excited  states  is  seen  to  reach  7 %  at  30,000  seconds. 
If  excited  state  energy  and  temperature  drops  are  taken  together,  one  can 
see  that  only  about  3%  of  the  initial  energy  has  been  radiated  by  the  end 
of  this  calculation. 
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Calculations  at  electron  densities  of  106  and  108  cm-3  show 
quite  different  behavior.  Both  of  these  cases  rapidly  reach  a  quasi¬ 
steady  configuration  (less  than  1000  seconds  required  at  106  cm-3)  and 
substantial  fractions  of  the  initial  energy  have  been  radiated.  Although 
the  curves  denoting  the  energy  fraction  in  excited  states  peak  a  bit  above 
10%,  they  decay  to  about  3.5%  by  the  time  one-half  the  initial  energy  has 
been  radiated  (-20,000  seconds  at  106  cm-3;  -3000  seconds  at  108 
cm-3).  In  these  cases,  the  infrared  (2P°-*2D°)  transition  is  responsible 
for  the  bulk  of  the  cooling.  Observe  that  in  both  Figures  14  and  15,  the 
2D°  state  density  peaks  at  about  13%  of  the  total  nitrogen  atom  density. 
This  common  feature  arises  because  the  radiative  transition  probability 
from  the  state  (see  Table  1)  is  quite  low,  so  collisional  excitation  and 
deexcitation  is  responsible  for  establishing  the  state  density  in  accord 
with  a  Boltzmann  distribution,  i.e.,  the  ratio  of  excited  state  densities 
in  Equation  13  becomes  independent  of  electron  density. 

Calculations  of  radiations  from  02  are  considered  next.  Once 
more,  the  02  density  has  been  equated  to  the  ion  and  electron  densities. 
Figures  16,  17,  and  18  contain  the  relevant  plots.  The  ground  and  excited 
states,  X3!’.  al&g*  and  k1!*,  have  been  abbreviated  as  3£,  ja,  and  in 
these  figures.  Also  note  that  the  present  calculations  considered  only 
two  optical  transitions,  as  mentioned  earlier. 

Figure  16  details  the  excitation  and  radiation  characteristics 
of  02  at  an  electron  density  of  104  cm-3.  From  the  first  parts  of  the 
figure,  it's  clear  that  the  state  requires  about  10,000  seconds  to 
attain  a  quasi-steady  state  but  the  ll  state  is  populated  virtually 
instantly  as  compared  to  the  time  scale  of  these  calculations.  The  radia¬ 
tion  rates  reflect  these  excitation  characteristics.  At  this  electron 
density,  much  less  than  a  percent  of  the  available  energy  goes  into 
excited  states,  and  only  about  2-1/2%  of  the  initial  energy  is  radiated 
during  the  calculation  time. 
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Figure  16.  Radiation  properties  of  Oj  at  an  electron  density  of  104  cm- 
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At  higher  electron  densities,  the  02  radiation  characteristics 
are  quite  different.  At  106  cm-3  (Figure  17),  the  quasi-steady  state  is 
attained  in  about  2500  seconds,  and  radiation  removes  nearly  one-half  the 
initial  energy,  as  indicated  by  the  temperature  curves.  At  most,  only 
about  5*  of  the  available  energy  is  tied  up  in  excited  states.  Note  that 
prior  to  -10,000  seconds,  radiation  from  the  1£+3£  transition  was  the 
dominant  cooling  mechanism,  but  afterwards  the  1a+3£  transition  dominates. 

Figure  18  illustrates  that  at  108  cm-3  electron  density,  the 
quasi-steady  state  is  obtained  instantaneously  as  compared  to  the  computa¬ 
tion  time  span.  The  1£-»-3£  transition  is  seen  to  be  dominant  over  the  full 
30,000  second  period  and  is  chiefly  responsible  for  the  energy  loss  over 
the  entire  computation.  About  1500  seconds  were  required  to  radiate  one- 
half  the  initial  energy.  Note  that  in  this  case,  radiation  cooling  lower¬ 
ed  both  ion  and  electron  temperatures  to  0.22  eV  after  30,000  seconds. 
This  particularly  low  temperature  could  be  achieved  because  the  metastable 
levels  of  02  are  quite  low  as  compared  to  the  metastable  states  of  0+,  0, 
N+,  and  N  (Refer  to  Figure  1).  Finally,  observe  that  in  this  calculation 
the  peak  fraction  of  available  energy  in  excited  states  was  about  7%. 

To  complete  the  discussion  of  these  three  cases,  it's  worth 
pointing  out  the  trend  as  to  which  of  the  02  transitions  dominates.  At 
lOf*  cm-3  electron  density  (Figure  16),  the  ll+3l  transition  is  dominant 
for  only  the  first  2000  seconds;  at  106  cm"3  (Figure  17)  it's  dominant  for 
the  first  10,000  seconds,  and  at  10®  cm-3  (Figure  18),  it  dominates  for 
the  entire  30,000  seconds.  This  trend  is  not  a  direct  temperature  effect 
(lower  ultimate  temperatures  at  higher  electron  densities)  because  the 
state  is  lower  energy  than  is  the  ll  state.  Instead,  the  trend  is 
explained  by  Equation  13,  i.e.,  it  is  the  play-off  between  collisional  and 
radiative  processes  which  is  important.  There  is  an  indirect  temperature 
effect  insofar  as  the  collisional  excitation  and  deexcitation  rate 
coefficients  are  temperature  dependent. 


The  discussion  of  single  species  systems  is  now  complete.  The 
remainder  of  this  section  will  be  devoted  to  calculations  based  on  species 
densities  and  initial  electron  and  ion  temperatures  extracted  from 
detailed  high  altitude  nuclear  burst  calculations.  In  these  calculations, 
the  time  scale  has  been  adjusted  to  span  the  most  interesting  phase  of  the 
radiation  cooling  process. 

The  Starfish  high  altitude  nuclear  test  event,  conducted  on  9 
July  1962  during  the  Fishbowl  series,  consisted  of  a  1.4  megaton  device 
detonated  at  400  km  altitude  above  Johnston  Island  in  the  Pacific. 
Recently,  Fajen25  has  performed  a  series  of  detailed  two-fluid  magneto¬ 
hydrodynamic  calculations  of  the  Starfish  event.  The  results  of  his 
computer  simulations  contain  the  best  available  spatially  and  temporally 
resolved  estimates  of  electron  and  species  densities,  ion  and  electron 
temperatures,  and  coupled  ion-neutral  dynamics.  Numerical  data  corre¬ 
sponding  to  90  seconds  after  detonation  have  been  extracted  from  Fajen 's 
calculations  and  have  been  used  as  the  time  equals  zero  initial  conditions 
for  radiation  cooling  calculations  which  are  presented  in  Figures  19  and 
20.  The  reader  may  wish  to  refer  to  Table  2  for  a  summary  of  these  data. 


Figure  19,  which  corresponds  to  an  altitude  of  395  km,  illus¬ 
trates  an  element  of  fireball  plasma  which  is  predominantly  atomic 
neutrals  (N  and  0)  at  a  total  atom  density  of  9.6  x  108  cm-3.  Oxygen  ions 
are  the  least  populous  species  (104  cm-3);  the  nitrogen  ion  density  is 
about  0.75%  of  the  neutral  atom  density.  Molecular  oxygen  has  a 
density  of  about  0.25%  of  the  neutral  atom  density.  The  plots  indicate 
that  with  the  exception  of  02 ,  the  initial  transients  required  to  populate 
the  excited  states  are  finished  in  a  few  tens  of  seconds.  As  one  might 
expect,  based  on  the  preponderance  of  oxygen  atoms  and  the  relatively  low 
temperature,  radiation  in  the  oxygen  red  line  (1D+3P  transition)  is  the 
major  cooling  mechanism;  the  electron  temperature  is  too  low  to  effec¬ 
tively  excite  the  oxygen  green  line  ( 1 S+ 1 D  transition).  The  second  most 
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Figure  19.  Radiation  properties  based  on  the  90  second  environment  at 
395  km  altitude  in  Fajen's  Starfish  calculation. 
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active  transition  (energy-wise)  is  the  infrared  2P°->-2D0  transition  in 
atomic  nitrogen,  but  after  about  50  seconds  (i.e.,  140  seconds  after 
detonation)  the  energy  loss  rate  is  about  a  factor  of  ten  below  that  for 
the  oxygen  red  line.  For  this  particular  calculation,  other  transitions 
and  species  are  relatively  unimportant  to  the  cooling  process. 

The  fraction  of  available  energy  tied  up  in  excited  states  is 
seen  in  Figure  19  to  peak  at  about  45*,  a  value  substantially  larger  than 
found  in  previous  calculations.  This  large  fraction  ocurs  because  the 
ratio  of  atoms  to  free  electrons  is  large  (-100),  and  the  electrons  have 
sufficient  thermal  energy  to  populate  the  excited  states.  Thus,  nearly 
half  the  avaialble  energy  is  rapidly  stored  as  excited  state  energy. 

The  calculation  depicted  in  Figure  20  corresponds  to  conditions 
extracted  from  Fajen's  Starfish  simulation  at  90  seconds  after  burst  at 
550  km  altitude.  In  this  case  the  electron  density  and  neutral  atom 
densities  are  nearly  equal.  The  0+  density  is  about  a  factor  of  2  larger 
than  the  N+  density,  while  the  0  density  is  about  a  factor  of  3  larger 
than  the  N  density.  Molecular  oxygen  is  a  minor  constituent.  Although 
both  ion  and  electron  temperatures  were  about  1.5  eV  initially,  very  rapid 
collisional  processes  caused  these  temperatures  almost  instantly  to  drop 
to  about  an  eV.  In  conjunction,  we  note  that  the  fraction  of  available 
energy  in  excited  states  peaks  at  nearly  18%  quite  early  in  the  calcula¬ 
tion.  It's  also  interesting  to  observe  that  except  for  N  and  02,  the 
population  densities  of  the  first  excited  states  of  the  species  reach  10% 
of  the  density  of  the  respective  ground  states.  For  N  and  02,  the  first 
excited  state  populations  reach  about  25%  of  the  respective  ground  state 
densities. 

The  combined  action  of  the  ionic  and  atomic  radiators,  all  of 
which  are  radiating  at  roughly  the  same  rate,  is  to  remove  about  59%  of 
the  initial  energy  over  the  1000  second  time  span.  (Keep  in  mind  that 
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Figure  20. 


Radiation  properties  based  on  the  90  second  environment  at 
550  km  altitude  in  Fajen's  Starfish  calculation. 
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these  calculations  do  not  include  effects  of  fluid  motion,  deionization 
chemistry,  etc.  which  would  act  in  a  complete  and  self-consistent  calcula¬ 
tion  during  the  1000  second  span  to  alter  the  available  thermal  energy 
and,  hence,  the  radiation  rate.)  Radiation  from  molecular  oxygen  is 
negligible,  from  the  cooling  standpoint,  because  the  O2  density  is  so  low. 

In  order  to  investigate  radiation  cooling  and  energy  partition 
long  after  a  nuclear  detonation,  we  consider  the  high  altitude  environment 
predicted  by  a  detailed,  two-fluid  numerical  simulation  26  of  a  megaton- 
range  explosion  at  200  km  altitude.  Calculations  depicted  in  Figure  21 
and  22  represent,  at  time  zero,  species  densities  and  temperatures 
extracted  from  the  detailed  simulation  at  60  minutes  after  detonation. 
These  radiation  computations  span  almost  3  hours. 

The  calculation  presented  in  Figure  21  represents  the  environ¬ 
ment  at  850  km  altitude  within  a  dense  portion  of  the  geomagnetic  field 
aligned  fireball /plume.  As  indicated  in  Table  2,  the  electron  density  is 
slightly  above  108  cm-3,  and  the  predominant  species  are  ions.  Ion  and 
electron  temperatures  were  both  initialized  to  0.504  eV  at  the  beginning 
of  the  radiation  cooling  calculation.  One  finds  that  this  system  almost 
instantaneously  attains  a  quasi-steady  state  with  only  a  bit  more  than  one 
percent  of  the  available  energy  associated  with  excited  states.  An  exam¬ 
ination  of  radiation  rate  curves  reveals  that  nitrogen  ions  are  the  most 
active  radiators.  In  fact,  N+  transitions  radiate,  over  10,000  seconds, 
more  energy  than  all  of  the  other  species  combined;  about  30%  of  the 
initial  energy  is  radiated  during  this  interval. 

The  final  case  to  be  considered  is  shown  in  Figure  22  and  repre¬ 
sents  conditions  extracted  from  the  aforementioned  simulation  at  60 
minutes  after  detonation.  This  set  of  species  and  temperatures  represents 
the  environment  at  1000  km  altitude  in  a  region  within  the  geomagnetic 
field  aligned  f ireball/plume  plasma  but  outside  of  its  dense  core.  Ionic 
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Figure  21 . 
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Radiation  properties  based  on  the  60  minute  environment  at 
850  km  altitude  in  a  calculation  of  a  megaton-range  nuclear 
explosion  at  200  km  altitude. 
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Figure  21  (Continued) 
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Figure  22.  Radiation  properties  based  on  the  60  minute  environment  at 
1000  km  altitude  in  a  calculation  of  a  megaton-range  nuclear 
explosion  at  200  km  altitude. 
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and  electron  temperatures  have  values  of  0.393  eV,  so  this  case  describes 
plasma  which  is  initially  much  cooler  than  that  in  previous  cases.  The 
radiation  cooling  curves  indicate  that,  at  these  low  temperatures,  nitro¬ 
gen  ions  and  oxygen  atoms  are  the  effective  radiators.  In  both  species, 
it  is  the  red  lines  which  account  for  the  cooling.  In  view  of  the  low 
temperatures  and  approximate  equality  between  electron  and  radiating 
species  densities,  it's  not  surprising  that  at  most  about  one  precent  of 
the  available  energy  was  contained  in  excited  states.  The  temperature 


curves  indicate  that  about  17%  of  the  initial  energy  was  radiated  over  the 
10,000  second  period. 


SECTION  4 

DISCUSSION  AND  CONCLUSIONS 

The  sample  calculations  of  the  preceding  section  rather  conclu¬ 
sively  demonstrate  that,  on  a  time  scale  of  hours,  forbidden  line  radia¬ 
tion  from  low  lying  metastable  states  of  atmospheric  constituents  can 
remove  substantial  fractions  of  the  available  thermal  energy.  From  the 
standpoint  of  computing  the  behavior  and  decay  of  fireball  plasma  from  a 
high  altitude  nuclear  explosion,  a  process,  such  as  radiation  cooling, 
which  significantly  changes  the  ion  and  electron  temperatures  may  be  quite 
important.  In  addition  to  altering  temperature  dependent  deionization 
chemistry  rate  coefficients  and  the  ion-electron  scale  height,  radiation 
cooling  may  result  in  a  reduction  in  the  diamagnetic  levitation  of  a 
plasma  in  the  geomagnetic  field.  This  last  effect  derives  from  the  fact 
that  ion-electron  thermal  energy  drives  a  system  of  electrical  currents 
within  the  plasma  which  interact  with  the  geomagnetic  field  and  cause  the 
plasma  to  move,  transverse  to  the  field,  toward  a  region  of  reduced  field 
( i . e . ,  generally  upward). 


The  results  presented  are  a  bit  artificial  in  the  sense  that 
they  represent  the  response  of  an  initially  cold  system  of  particles  to  a 
delta-function  heating  pulse  at  zero  time.  In  a  more  realistic  calcula¬ 
tion,  the  system  at  some  non-zero  initial  temperature  would  be  subject  to 
non- instantaneous  heating  (and  cooling)  as  the  result  of  hydrodynamic 
evolution.  Nevertheless,  the  results  provide  a  basis  for  estimating 
e-folding  time  constants  for  cooling  and  for  particular  forbidden  line 
radiation  rates  under  a  variety  of  conditions.  Although  most  of  these 
forbidden  lines  are  in  the  visible,  a  few  are  in  the  ultraviolet  and  short 
wavelength  infrared. 


Earlier,  in  Section  2,  it  was  pointed  out  that  excited  state 
populations  and  radiation  rates  do  not  scale  in  a  simple  way  with  electron 
density.  The  mathematical  statement  of  the  scaling  under  quasi-steady 
conditions  was  presented  in  Equation  13,  and  the  computed  results  nicely 
demonstrate  the  variation  of  excited  populations  with  electron  density 
under  more  general  conditions.  The  intertwining  of  temperature  dependent 
excitation  and  deexcitation  rate  coefficients  and  radiative  transition 
probabilities  makes  it  difficult  to  sort  out  the  significance  of  individ¬ 
ual  line  radiations  relative  to  the  others  without  doing  quite  a  bit  of 
work. 


A  very  interesting  result  of  the  calculations  concerns  the 
partition  of  available  energy  between  excited  state  and  thermal  energies. 
The  computations  indicate  that  when  the  temperature  is  an  electron  volt  or 
lower,  and  when  the  electron  density  is  comparable  to  or  exceeds  the 
density  of  radiating  species,  then  the  fraction  of  available  energy  tied 
up  in  excited  states  is  -10%  or  less.  However,  if  the  radiating  species 
density  substantially  exceeds  the  electron  density  (e.g.,  far  more  N  or  0 
atoms  than  electrons),  then  a  large  fraction  of  the  available  energy  may 
go  to  producing  excited  states,  even  if  the  temperature  is  fairly  low 
(<0.4  eV).  Under  such  circumstances,  fluid  dynamical  calculations  which 
compute  ion-electron  pressures  without  accounting  for  thermal  energy  which 
would  be  absorbed  in  the  production  of  excited  states  will  possibly 
miscompute  pressure  gradient  forces  and  diamagnetic  effects  by  perhaps  a 
factor  of  two.  One  remedy  for  this  problem  would  be  to  carry  explicitly 
excited  state  populations.  However,  the  present  calculations  indicate 
that  at  electron  densities  of  greatest  interest  (106  cm-3  and  larger),  the 
quasi-steady  state  is  attained  rather  rapidly.  Therefore,  it  should  be 
possible,  in  all  but  highly  dynamic  situations,  to  compute  by  the 
quasi-steady  approximation  the  excited  state  populations  without  explic¬ 
itly  carrying  them.  Thus,  standard  fluid  treatments  could  be  improved  by 
computing  the  energy  tied  up  in  excited  states  and  by  subtracting  this 
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energy  from  the  thermal  energy  prior  to  a  calculation  of  ion  electron 
pressure  or  temperature. 

In  summary,  it  has  been  determined  that  high  altitude  nuclear 
burst  simulations  which  span  substantial  time  intervals  (hours)  need  to 
include  radiation  cooling  by  forbidden  line  radiation.  The  quasi-steady 
state  approximation,  which  avoids  the  necessity  of  carrying  excited  state 
densities  explicitly,  appears  to  be  acceptable  for  all  but  highly  dynamic 
problems.  Furthermore,  it  has  been  shown  that  failure  to  account  for  the 
partition  between  excited  state  and  thermal  energies  can  introduce 
non-trivial  errors  into  pressure  and  other  temperature  dependent  terms  of 
the  magnetized  fluid  equations. 
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Transition 


Deexcitation  Rate  Coefficient 
(cm3 /sec) 


1.3  x  10-8  (Tg)-1/2 

6.3  x  10-9  (Te)"1/2 

2.4  x  10-8  (Te)-1/2 


t  From  reference  13, 


TABLE  A-2.t  DEEXCITATION  RATE  COEFFICIENTS  FOR  THE  LOW-LYING  METASTABLE 
STATES  OF  N4  (Te  in  eV) 


Transition 

Deexcitation  Rate  Coefficient 
(cm3 /sec) 

lD  -  3P 

4.8  x  10-8  (Tg)-1/2 

lS  -  3P 

3.2  x  10-8  (Tg)"1/2 

lS  -  lD 

3.3  x  IQ-8  (Te)-1/2 

*  From  reference  13 


f 


TABLE  A-3. 


Te(eV) 


ELECTRON  IMPACT  EXCITATION  AND  DEEXCITATION  RATE  COEFFICIENTS 
FOR  THE  LOU-LYING  STATES  OF  ATOMIC  OXYGEN. 


1.92  (-18)“ 
1.1  (-  9)b 

5.28  (-14) 
1.7  (-  9) 


1.76  (-12) 
2.2  (-  9) 

3.28  (-11) 


1.78  (-28) 
2.1  (-  9) 

1.96  (-19) 
2.1  (-  9) 

2.10  (-16) 
2.1  (-  9) 

6.04  (-14) 
2.3  (-  9) 


2.25  (-19) 

2.2  (-13) 

1.76  (-14) 
3.9  (-11) 

8.07  (-13) 

2.3  (-10) 

1.52  (-11) 
8.7  (-10) 


0.7 

1.21 

(-101 

1  7.25 

(-13] 

1  5.4  ( 

3.5 

1-  9] 

1  2.5 

I-  9| 

1  1.5  ( 

1.0 

3.43 

(-1°| 

1  4.93 

( -3.2 1 

1  1.38  ( 

4.4 

(-  9] 

1  2.9 

(-  9j 

1  2.3  ( 

5.20  (-10) 

4.8  (-  9) 

7.94  (-10) 
5.2  (-  9) 

1.21  (-  9) 

5.8  (-  9) 

1.84  (-  9) 


1.06  (-11) 

3.1  (-  9) 

2.32  (-11) 
3.4  (-  9) 

5.15  (-11) 
3.7  (-  9) 

1.16  (-10) 

4.2  (-  9) 


(-  9) 

2.21 

(-10) 

(-  9) 

4.6 

(-  9) 

(-  9) 

2.8 

(-10) 

(-  9) 

4.6 

(-  9) 

(-  9) 

3.3 

(-10) 

(-  9) 

4.5 

(-  9) 

2.58  (-  9) 

5.3  (-  9) 

2.23  (-  9) 

4.4  (-  9) 


3.8  (-10) 
4.5  (-  9) 

3.7  (-10) 
4.1  (-  9) 


1.97  (-10) 

2.7  (-  9) 

2.76  (-10) 

3.1  (-  9) 

3.98  (-10) 

3.6  (-  9) 

5.30  (-10) 

3.9  (-  9) 

7.16  (-10) 

4.5  (-  9) 

7.76  (-10) 

4.6  (-  9) 

7.5  (-10) 

4.2  (-  9) 

6.9  (-10) 

3.7  (-  9) 

6.4  (-10) 

3.4  (-  9) 


Numbers  In  parenthesis  indicate  the  power  of  ten  by  which  the 
enterles  are  multiplied. 

Excitation  rate  coefficient. 

The  corresponding  deexcitation  rate  coefficient. 

From  reference  13. 
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TABLE  A-4.  ELECTRON  IMPACT  EXCITATION  AND  DEEXCITATION  RATE  COEFFICIENTS 
FOR  THE  LOW-LYING  STATES  OF  ATOMIC  NITROGEN. 


Numbers*  p  enthesls  indicate  the  power  of  ten  by  which  the 
enterles  «.i ~  multiplied. 

Excitation  rate  coefficient. 

The  corresponding  deexcitation  rate  coefficient. 

From  reference  13. 


Numbers  in  parenthesis  indicate  the  power  of  ten  by  which  the 
entries  are  multiplied. 
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